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ABSTRACT 

The implications of hyperon pairing for the thermal evolution of neutron stars con- 
taining hyperons are investigated. The outcome of cooling simulations are compared for 
neutron star models composed only of nucleons and leptons, models including hyperons, 
and models including pairing of hyperons. We show that lambda and neutron pairing 
suppresses all possible fast neutrino emission processes in not too massive neutron stars. 
The inclusion of lambda pairing yields better agreement with X-ray observations of pul- 
sars. Particularly, the surface temperatures deduced from X-ray observations within 
the hydrogen atmosphere model are more consistent with the thermal history of neu- 
tron stars containing hyperons, if the critical temperature for the onset of lambda and 
nucleon pairing is not too small. 
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1. Introduction 

Relativistic calculations of the composition of neu- 
tron star matter lead to the conclusion that neutron 
stars are composed not only of nucleons and lep- 
tons but also of hyperons and, possibly, of nucleon 
isobars (see, e.g., Pandharipande 1971, Glendcnning 
1985, Schaffner & Mishustin 1996, Balberg & Gal 
1997, Huber et al. 1997). It was proposed by Prakash 
et al. (1992) that the presence of hyperons might lead 
to rapid cooling even if the proton fraction is too small 
to allow for the nucleon direct Urea process and no ex- 
otic states (like meson condensates and quark gluon 
plasma) exist (see also Prakash 1994). This state- 
ment was confirmed by detailed cooling simulations 
(Haensel & Gnedin 1994, Schaab et al. 1996). The 
role of hyperons in the cooling history of a neutron 
star is even enhanced by the generally accepted neu- 
tron pairing in the interior of the star (see, e.g., Takat- 
suka 1972, Amundsen & 0stgaard 1985). Whereas 
the direct nucleon Urea is suppressed by neutron pair- 
ing (Lattimer et al. 1994, Page 1995) most of the 
hyperon processes are not. Haensel & Gnedin (1994) 
and Schaab et al. (1996) found that unsuppressed 
rapid cooling due to hyperon induced neutrino emis- 
sion processes results in surface temperatures that are 
too low compared with soft X-ray and extreme UV 
observations. 

Recently Balberg & Barnea (1998) calculated the 
^So gap energy of A hyperons in neutron star matter. 
The obtained gap energy is similar to the "'^Sq gap en- 
ergy of protons. As we will show here hyperon pair- 
ing considerably changes the cooling behavior of a 
neutron star containing hyperons (which we will call 
hyperon stars in the following) and leads to better 
agreement with observations. 

In this LeMer we study the effect of hyperon pairing 
on the thermal evolution of hyperon stars and com- 
pare the results with "conventional" neutron stars, 
which are composed of nucleons and leptons only, 
and with observations. We begin with the underlying 
physics by discussing the relevant neutrino emission 
processes and the pairing of hyperons (Sect. 2 and 
3) and describing in Sect. 4 the relativistic equations 
of state (EOS's) used. The results of the simulations 
are discussed in Sect. 5. In Sect. 6 we compare our 
results with observations and end with some conclu- 
sions. 



2. Neutrino Emission Processes 

The early evolution of hot neutron stars is com- 
pletely dominated by the cooling via neutrino emis- 
sion. Only after about 10'"' yrs does the photon radi- 
ation from the star's surface take over and dominate 
the late evolution. This late cooling might be slowed 
down by various heating processes (see Schaab et al. 
1998b and references therein). The neutrino emis- 
sion processes can be divided into slow and enhanced 
processes depending on whether one or two baryons 
participate. Due to the rather different phase spaces 
associated with both kind of processes the emission 
rates differ by several orders of magnitude. 

In hyperon stars various enhanced neutrino emis- 
sion processes can occur, as listed in Tab. 1. The 
processes are only noted in one direction. Because 
of /3-equilibrium the inverse reaction (with 9 replaced 
by v) occiirs at the same rate. The first process is 
the direct nucleon Urea process, the remaining pro- 
cesses are comprised by the term "direct hyperon 
Urea". The relative emissivities TIbiB2 = ^v{Bi 
B2 + e~ + D)/ev{n ^ p -|- e^ -|- i^) of these processes 
are given with respect to the emissivity of the direct 
nucleon Urea process (Prakash et al. 1992). Simulta- 
neous conservation of energy and momentum requires 
that the triangle inequality, < +Pe > ^'^^ the 
two inequalities obtained by cyclic permutation are 
fulfilled for the Fermi momenta pf. If the inequalities 
are not fulfilled the process is extremely unlikely to 
occur and the corresponding emissivity vanishes. The 
availability of the various fast processes is thus de- 
pendent on the partial concentrations of each baryon 
species, which in turn is determined by the EOS. 

If one of the participating baryons pairs in a su- 
perfluid state the emissivity is suppressed by an ap- 
proximately exponential factor: T^gf ~ cxp(— cTc/T) 
for T <C Tc where is the critical temperature and c 
is a constant of the order unity. Since the Fermi mo- 
menta of the hyperons are sufficiently small they pair 
in a ^ So state (see Sect. 3), as is the case for protons at 
slightly lower densities. Due to their high Fermi mo- 
menta, neutrons are expected to pair in an anisotropic 
^P2 state. We refer to Levenfish & Yakovlev (1994) 
for fitting formulas of T^sf, that are valid over the 
whole temperature range T < Tq for both types of 
pairing state. 

Table 1 summarizes all possible direct Urea pro- 
cesses with the corresponding relative emissivities. 
Also included are the critical density ric and the 
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Table 1: Possible direct Urea processes in hyperon stars 

TMl-ml 



RHF8 





Process 




T^-BiSa "factor 








Mth/Mo 
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^ p + e" 


+ v 


1 


1.31 


0.78 


1.96 


1.09 


E- 


^ A + e" 


- +v 


0.2055 


1.89 


1.30 


2.99 


1.35 


A 


— >■ p + e~ 


+ V 


0.0394 


1.90 


1.31 


3.11 


1.38 


E- 


■ — » n + 


+ i> 


0.0125 


2.37 


1.45 


2.44 


1.22 




A + e" 


' +P 


0.0175 


2.65 


1.52 


6.80 


1.65 


E- 


^EO+e 


- +iy 


0.6052 


4.03 


1.60 


7.43 






^ EO + e- 


~ +i> 


0.0282 


4.03 


1.60 








^ SO + e- 


- +v 


0.2218 


4.27 


1.61 






S° ■ 


E+ + e- 


~ +v 


0.0564 


4.30 









Note. — "R-BiB^ is the relative cmissivity ev{B\ B2 +e +v)/tv{r\ — > p + e Uc the critical density above which the process 

is allowed, and M^jj the corresponding threshold mass of a non-rotating neutron star (both values depend on the underlying EOS, here: 
TMl-ml and RHF8). 



threshold mass of a non-rotating neutron star Mth 
above which the respective process is allowed for two 
representative EOS's (see Sect. 4). Besides the direct 
nucleon Urea process n ^ p-)-e^ -l-P, the most impor- 
tant processes are E~ ^ A -I- e~ -I- z?, A ^ p -I- + 
E~ n + e~ + V, and K + e~ + v. Other 

processes are only possible at very high densities just 
below or above the maximum central density of a 
stable hyperon star. As a result, the gap energy 
of the A particle is the most important ingredient 
(following those of the nucleons) for simulating the 
cooling history of the star, whereas the gap ener- 
gies of the other hyperons need to be known only 
for the most massive stars. The corresponding mod- 
ified hyperon Urea processes obtained by adding a 
bystander baryon are only important if no enhanced 
process is allowed since their emissivities are consid- 
erably smaller (see Maxwell 1987). 

3. Superfluidity of Hyperons 

The presence of an attractive two-body A — A in- 
teraction is implied by doubly-strange hypernuclei, 
where two A hyperons are trapped in a single nucleus 
(Imai 1992). The properties of the decay products of 
such hypernuclei indicate that the separation energy 
of the two A hyperons, -Baa, is larger than twice the 
separation energy of a single A from the same core 
nucleus, 2Ba- The bond energy derived from these 
experiments is A_Baa = ^aa — 2-Ba ~ 4 — 5 MeV, 
somewhat less than the corresponding bond energy 



of nucleons in nuclei, Bmn ~ 6 — 7 MeV. Hence, A 
hyperon pairing similar to that of nucleons is indeed 
expected in high density matter. 

In our simulations presented below we use the A 
hyperon ^Sq gaps as estimated recently by Balberg & 
Barnea (1998). Their AA potential is based on the 
G-Matrix parameterization (Lanskoy & Yamamoto 
1997), which accurately reproduces the bond energy 
of various doubly strange hypernuclei. In this model, 
the pairing energy is calculated as a function of the A 
Fermi momenta and the background baryon density, 
while neglecting relativistic and polarization effects. 

In general, A pairing was found to extend from the 
threshold density for A formation to the density where 
the A Fermi momenta exceeds 1.3 fm~^ (typically cor- 
responding to a A fraction of 0.15 — 0.2). The corre- 
sponding gap energies were found to be several tenths 
of a MeV, where the exact value depends on both the 
A fraction and the background matter density. These 
values are actually similar to those found for proton 
^So pairing with a neutron matter background (El- 
gar0y et al. 1996), again indicating the general sim- 
ilarity between hyperon-related and nucleon related 
interactions expected at supernuclear densities. 

The pairing gaps are also strongly dependent on 
the effective mass of the particles on the Fermi sur- 
face, since the pairing interaction competes with the 
particles' kinetic energy. Relativistic calculations typ- 
ically predict that at supernuclear densities the the 
baryon effective masses are considerably lower than 
the bare masses, but the exact value on the Fermi 
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Fig. 1 — Critical temperatures for the TMl-ml EOS. 



surface is highly model dependent. We thus treat the 
value of the effective mass as a parameter, which al- 
lows to examine the sensitivity of our results to the 
uncertainty in the A gap energies. 

Figure 1 shows the critical temperature for A 
"'^So pairing expected for the equilibrium composi- 
tions found with the TMl-ml-model EOS (Schaffner 
& Mishustin 1996). The critical temperature is re- 
lated to the gap energy A(0) at zero temperature by 
Tc = aA(0)/A;B, where the proportional factor a de- 
pends on the pairing type (see Levenfish & Yakovlev 
1994). Note the strong dependence of the critical tem- 
perature on the A effective mass, m^, (unlike the 
density range in which ^Sq paring exists, which is 
practically independent of m^). We show the criti- 
cal temperatures when the A mass on the Fermi sur- 
face is taken as equal to the bare one, m\/m^ = 1, 
and also when this mass is reduced by a factor of 
m\/m\ = 0.7. The latter case is a reasonable repre- 
sentation for the value of the bulk effective mass as 
found in the sclfconsistent calculation of the TMl-ml 
model. We stress again that the mass on the Fermi 
surface may differ considerably from the bulk effec- 
tive mass, and the two cases will be used here mainly 
to examine the dependence of the thermal history on 
the A critical temperature. Also shown are the critical 
temperatures for proton -'^So and neutron pairing, 
as calculated by Wambach et al. (1991) and Amund- 
sen & Ostgaard (1985), respectively. Note that neu- 
tron pairing exists over the entire supernuclear 
density range in this specific model. 



4. Equations of State 

The EOS is quite unknown at high density. For 

the present study, we choose two recent rclativistic 
parameterizations, which have been proven to be able 
to describe properties of nuclei and/or nucleon scat- 
tering data. 

The first set (TMl-ml) is a rclativistic mean-field 
calculation (model 1 of Schaffner & Mishustin 1996). 
The nucleon parameters arc adjusted to properties of 
nuclei (set TMl, see Sugahara & Toki 1994). The 
vector potential is nonlinear in density to get a soft 
behavior at high density in accordance with Dirac- 
Briickner calculations. The hyperon coupling con- 
stant are fixed by SU(6) symmetry and hypernuclear 
data. 

The other set (RHF8 from Huber et al. 1997) 
uses Relativistic-Briickner-Hartree-Fock results up to 
2-3 times normal nuclear density. Hypcrons arc im- 
plemented within the Relativistic-Hartree-Fock ap- 
proach. The EOS calculated within Relativistic- 
Hartree-Fock are usually much softer than those de- 
rived from rclativistic mean-field models due to the 
additional degrees of freedom. Again, SU(6) symme- 
try as well as hypernuclear data are used to fix the 
hyperon coupling constants. 

The hyperons appear at a slightly higher density 
for RHF8 than in TMl-ml except for the S" which 
appears around the same density of 2po- The maxi- 
mum fraction of S" is higher in RHF8 than in TMl- 
ml, while the A fraction is about a factor of two 
smaller, This has an impact on the onset of the vari- 
ous hyperon cooling processes as seen in Tab. 1. The 
critical density for hyperon induced cooling processes 
is usually lower for TMl-ml than for RHF8, except 
for the process S~ n + e~ + v, which is however 
suppressed by neutron pairing. 

5. Results 

The general rclativistic equations of stellar struc- 
ture and thermal evolution (see Thorne 1977) were 
numerically solved via an implicit finite difference 
scheme by a Newton-Raphson algorithm (see Schaab 
et al. 1996 for more details)^. The surface tempera- 
tures deductxl from the observed thermal X-ray spec- 
Tables with detailed references to the used ingredients, the used 
observational data, and the obtained cooling tracks can be 
found on the Web: http://www.physik.uni-muenchen.de/sek- 
tion/suessmann/astro/cool/schaab.0897. 
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Fig. 2. — Cooling of neutron and hyperon stars con- 
structed for the TMl-ml EOS. The observational 
data obtained within the magnetic hydrogen atmo- 
sphere model are marked by hollow circles, whereas 
the one obtained within the pure iron atmosphere 
model are marked by solid ones. 



tra depend on the used atmosphere model. We plot 
the surface temperatures along with their 2a error 
bars for the two extreme cases of a pure iron atmo- 
sphere (solid circles in Figs. 2 and 3) and of a mag- 
netic hydrogen atmosphere (hollow circles, see Schaab 
et al. 1998a for more details). The deduced surface 
temperatures can only be understood as upper limits 
to the evolution temperature, since magnetospheric 
emission may irradiate the star's surface. On the 
other hand, the magnetic hydrogen atmosphere model 
is likely to underestimate the actual surface temper- 
ature (Potekhin ct al. 1997). The two extreme val- 
ues of the respective star's surface temperature should 
therefore give a good estimate of the maximum range 
of the true evolution temperature. 

Fig. 2 shows the cooling behavior of hyperon stars 
constructed for the TMl-ml EOS. The surface tem- 
perature as measured by a distant observer is plot- 
ted against the star's age for different star masses 
(a larger mass corresponds to a larger central density, 
thus crossing the thresholds of more cooling processes; 
see Lattimer et al. 1991). 



• M — 0.7Mq: This model is an example for a 
slow cooling neutron star. The interior does 
not contain hyperons, nor is the proton frac- 
tion sufficiently high to allow for the (lirc;c;t nu- 
cleon Urea. Note that the relatively low surface 
temperature is caused by the inclusion of the 
superfluid pair breaking and formation process 
(Voskresenskii & Senatorov 1987, Schaab et al. 
1997), which is not included in former works 
(e.g. Umeda et al. 1994, Page 1995, Schaab 
et al. 1996). 

• M = 1.2Mq: The interior consists still only of 
nucleons, but the nucleon direct Urea is now 
possible. This yields a faster cooling which is 
slowed down by neutron pairing below the crit- 
ical temperature Tc w 8 x 10^ K. 

• M = IAMq: At n w 1.8 - 1.9no S" and 
A hyperons begin to populate the central re- 
gion of the star. Besides the direct nucleon 
Urea the processes E" A -|- e^ -|- P and 
A p + e~ + 9 are also possible. The two 
direct hyperon processes are suppressed by A 
pairing below the critical temperature T^. At 
the outer boundary of the hyperon core Tc is 
rather small (« 10® K for m,\/m\ = 1 and 
Ri 10^ K for m\/mA = 0.7). The hyperon pro- 
cesses are therefore significantly suppressed only 
in the innermost part of the star. The surface 
temperature of a middle aged, enhanced cooling 
neutron star (lO^yrs ^ t ^ lO^yrs) is directly 
related to the critical temperature of its super- 
fluids (Lattimer et al. 1994, Page 1995, Schaab 
et al. 1996) and depends therefore on the con- 
sidered gap energy model. For the A gap model 
with mJ^/mA = 1, the critical temperature of 
A pairing is larger than for neutron '^P2 pair- 
ing. In this specific model, the cooling behavior 
is therefore determined by the critical tempera- 
ture for the latter. 

• M = I.6M0: At densities n ^ 3no the scalar 

AA interaction becomes repulsive in our model 
TMl-ml, and the A's no longer pair in a 
^So state. The direct hyperon Urea processes 
are therefore unsuppressed and the interior of 
the star cools very rapidly. This result relies 
however on the specific assumptions that we ne- 
glected pairing of other hyperons (S°'^, 'EP'~) 
and pairing of A in a '^P2 state. If such further 
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Fig. 3.— Same as Fig. 2 but for the RHF8 EOS. 

pairing is allowed the cooling behavior should 
be similar to that of a M = IAMq star. 

Fig. 3 shows the cooling behavior for the RHF8 
EOS. As for TMl-ml we have four scenarios corre- 
sponding to the four star masses. The cooling be- 
havior of these scenarios are quite similar for both 
cases of EOS's, although some aspects arc somewhat 
model dependent (TMl-ml stars generally cool more 
rapidly). We emphasize the effect of A pairing on the 
thermal evolution for both EOS's, which is explicit in 
the difference between the M = IAMq (M = 1.5Mq 
for RHF8) thermal sequences and the M = 1.6Mq se- 
quences, as the latter essentially assume no hypcron 
pairing. This effect is enhanced, of course, for larger 
A pairing energies. 

6. Conclusions and Discussion 

In this Letter we have studied the effect of hy- 
peron pairing on the thermal evolution of hyperon 
stars. Our general finding is that A-hyperon pairing 
along with nucleon pairing is sufficient to suppress 
fast hyperon-induced cooling processes in all but the 
most massive neutron stars, in a similar fashion as 
the suppression of fast nucleon processes by nucleon 
pairing. Pairing of other hyperons becomes impor- 
tant only for cooling of the most massive stars with 



M > 1.6 Mq. 

As a result of this suppression, the surface tem- 
peratures obtained in cooling simulations of hyperon 
stars are considerably higher than those fomid for 
the same stars when the A hyperons are taken to 
be in a normal state. In particular, this allows bet- 
ter compatibility with surface temperatures deduced 
within the magnetic hydrogen atmosphere model from 
observed thermal emission from pulsars, with which 
cooling of hyperon stars through unsuppressed rapid 
cooling has been previously found to be inconsistent 
(Haensel & Gnedin 1994, Schaab et al. 1996). 

Our quantitative results demonstrate a significant 
sensitivity of the surface temperatures to the mod- 
eling of the hyperon pairing. We also call attention 
to the dependence of the surface temperature on the 
details of the high density EOS through the equilib- 
rium concentrations of the baryon species. In this 
sense, comparison of the outcome of theoretical cool- 
ing models with soft X-ray and extreme UV obser- 
vations provides a powerful tool for investigating the 
interiors of pulsars. Future measurement of pulsar 
thermal emission may indeed offer valuable indica- 
tion regarding the physics and composition of very- 
high density matter, although much improvement is 
still necessary in both theory and observation. 

Our main conclusion is that that the existence of 
hyperons in neutron stars cannot be excluded by com- 
paring their predicted thermal history with actual ob- 
servations from pulsars, provided the hyperons are al- 
lowed to form a superfluid state - similar to that of the 
nucleons. Though we have neglected several effects, 
such as medium effects on the modified Urea process, 
possible meson condensation and accreted envelopes, 
we believe this main conclusion to be quite robust. 
The study of these effects, as well as pairing of other 
hyperon species and other hyperon coupling models 
will be addressed in future work. 

Ch. S. gratefully acknowledges the Bavarian State 
for financial support. J. S.B. is supported by the 
Alexander von Humboldt-Stiftung and by the U.S. 
Department of Energy under Contract No. DE-AC03- 
76SF00098. 
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